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A new analytical method of performing ERBL evolution is described. The main goal 
is to develop an approach that works for distribution amplitudes that do not vanish at 
the end points, for which the standard method of expansion in Gegenbauer polynomials 
is inefficient. Two cases of the initial DA are considered: a purely flat DA, given by 
the same constant for all x, and an antisymmetric DA given by opposite constants for 
X ^ 1/2. For a purely fiat DA, the evolution is governed by an overall (xx)* dependence 
on the evolution parameter t times a factor that was calculated as an expansion in t. For 
an antisymmetric flat DA, an extra overall factor |1 — 2a;p' appears due to a jump at 
X = 1/2. A good convergence was observed in the t < 1/2 region. For larger t, one can 
use the standard method of the Gegenbauer expansion. 

Keywords: ERBL evolution; flat distribution amplitude; analytic methods. 

PACS numbers: 12.38.-t, ll.lO.-z 

1. Introduction 

Recent BaBar datcPon the 7*7 —> 7r° transition form factor correspond to approx- 
imately logarithmic InQ^ raise of the combination Q^F^»^^o{Q'^) in the region of 
very high momentum transfers 10 to 40 GeV^, where perturbative QCD approaclP 
predicts nearly constant behavior for this combination. It was proposed^ to ex- 
plain the BaBar "puzzle" by assuming that the pion distribution amplitude (DA) 
is "flat": ipTr{x) — fn- 

In general, DAs ipTr{x,iJ,) depend on the normalization scale /z, tending to the 
"asymptotic" shape y:>^{x) = 6fT^x{l — x) in the /it — >■ 00 limit, for any intial pion 
DA. A standard wajEEIto this result is to expand the initial pion DA over the eigen- 
functions a;(l — x)Cn^^{'2x — 1) of the evolution kernel. Each Gegenbauer projection 
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then changes as [lnln(/io/A)/ lnln(/i/A)]"''"/^'' when /i increases. All anomalous di- 
mensions 7„ are positive, except for 70 which is zero, hence only the ~ — x) 
part survives. For a pion DA given by a sum of a few Gegenbauer polynomials, 
this method gives a convenient analytic expression for the DA evolution. However, 
to get a flat pion DA, one should formally take an infinite number of Gegenbauer 
polynomials (this difficulty was mentioned in a recent papeii^ , and one should take 
a very large number of terms to get a reasonably precise (point by point) result for 
the evolved DA. In fact, a long time ago, it was establishecPthat, if the initial pion 
DA ipnix, /io) has a ~ [x(l — x)Y form near the end points a; = 0, 1, with a small r, 
then evolution just increases this power: 

r ^ rip) - r(Mo) + 2CFHH^i/ K) / Mp^/ A)]/ . 

This result was obtained by incorporating the large-n behavior 7„ ~ AC ph\n of 
the anomalous dimensions in QCD, and summation of the leading end-point region 
terms. Thus, one may expect that the evolved form of the flat pion DA should be 
close to [x{l-x)Y'^'' in[in(M/A)/in(/'o/A)]//3o^ particular, this statement was repeated 
in a more recent papei^ on the basis of an analysis similar to that of Ref. [5] . 

Our goal is to develop a method that allows one to get explicit analytic expres- 
sions for the evolution of the flat distribution amplitude in QCD. As we will see, our 
results not only confirm the above expectation, but also provide a regular way of 
calculating corrections to it. We also apply this new method to obtain the evolution 
of the DA d{x) ^ sgn(a; — 1/2) that is antisymmetric with respect to the central 
point X — 1/2. Such DA may correspond to the so-called "D-term"'Sl in generalized 
parton distributions. 



2. Evolution of Flat DA 
2.1. Evolution equation 

Evolution of the pion DA in the leading logarithm approximation is governed by 

dip{x, t) 



dt 



[y{x,y)]+ip{y,t)dy, 



(1) 



where t — 2Cf lnln(/x/A)/feo is the leading logarithm QCD evolution parameter, 
and \V{x, y)\+ is the evolution kernel which has the property of a plus distribution 



[y{x,y)]+ = V{x,y) - 5{y - x) / V{z,y)dz 



with respect to its first argument. In our case, 

1 



V{x,y)^ 



^(1 

y 



y — X 



e{x <y) + 



^(1 



1 



x-y 



0{y < x) , 



(2) 



(3) 



where we use a standard notation x = \ — x and y ^ I — y. One may rewrite the 
evolution equation as 

d(p{x,t) '■^ 



dt 



[V{x, y)ip{y, t) - V{y, x) ip{x, t)] dy 



(4) 
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It is clear that the singularities of V{x,y) and V{y,x) at x = y cancel each other 
in the integrand above, even though the subtraction there does not look like a 
"+" -prescription with respect to the integration variable. Adding and subtracting 
V{x, y)(p{x) in the integrand, we obtain the equation 



d^ix,t) 
dt 



I V{x,y)[(p{y,t) - (p{x,t)]dy - (p{x,t) [V {y , x) - V {x , y)]dy 
Jo 



(5) 



in which the first term has the structure of the "+" -prescription with respect to the 
integration variable, so that l/{x — y) singularity of V{x,y) is cancelled by zero of 
(fiiyjt) — (p{x,t) aX X — y. The integral in the second term (call it —v{x)) is also 
finite. Taking the Ansatz 



99(,x,t)=e*''(-)$(x,<) , 
we obtain for the equation 



V{x,y)[e'i^^y^--^-U{y,t) - Hx,t)]dy, 



(6) 



(7) 



which does not have the second term. The solution for ^{x,t) may be written as a 
series in t 



n=0 



with the functions $„(x) satisfying the recurrence relation 



^n+i{x) = [ V{x,y) 
Jo 



.1=0 ^ ' 



(8) 



dy. (9) 



2.2. Singular Part 

It is instructive to consider first an auxiliary situation when the evolution kernel is 
given by the singular part only 



V^^^{x,y) 



My-x) 

of the QCD kernel. In this case 



0{x < y) + (x x, y^y) 



^;«ing(a;) = _ /" [y''i"g(y, x) - ^^'"^(a;, y)]dy = 2 + In(xx) 
Jo 

and the recurrence relation is given by 

^„.r(.) = [ vix, y) \± j^^iiy) (m I)""' - 



(10) 



(11) 



dy , (12) 
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(13) 

{x)= j y(x,y) [$o(y)ln(^) +<fi(y)-<i>i(x)l dy, (14) 

In L \XX/ 1 

$2(2/) - $2 (x)] dy.{lb) 



yy_ 

XX 



or, explicitly for the first terras, 

$i(a;) = / V{x,y)[<i>o{v)~Mx)]dy, 
Jo 
fi 

$2(2;) = 

/o 

$3(2;)= / y(a;,y) [$o(y)ln'(^) +2<i>i(2/)ln( 
Jq l \xx/ \ 

In this approximation, we can write 
If we take the flat DA for t — 0, i.e., $(x, 0) = 1, this gives 

Mx) = 1 

$i(x) = 

^2{x) = -21nxlnx 

<i>3 (x) — 3 ln(a;a;) In a; In a; + 2 [In a; Li2 (a;) + In a; Li2 (a;)] 
-4 [Li3(a;)+ Li3(5)]+8C(3) 
The graphical results for the expansion components are given in Fig.l 



(16) 

(17) 
(18) 
(19) 

(20) 









0.5 



(a) (b) 
Fig. 1. Expansion components <I>2'"^(^') {^) ^^'^ ^s'^^i^) (b)- 

As far as J^[V{x,y)]-^-dx = 0, evolution does not change the normalization 
integral for ip{x, t). In particular, if we expand f{x, t) in t 



— ^ ri ' 



n=0 



we should have 



(Pn{x) dx = SnQ 



(21) 



(22) 
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However, in our Ansatz ip{x,t) = e*'-^"'"'"^^^^)$(a;, i), the (tlnxx)^ terms are 
summed to ah orders, while the series over fnix) is restricted to some finite order 
N. As a result, the approximants (/3(Ar)(a;, t) are not normalized to 1. In particular, 
if we keep the terms up to $2(2;), the normalization integral is given by 

hit) ^ £ ^^^f{x, t) dx = e^* Y§T^{ ^ ^^+^*)' " ^^(^ + }^^^) 

with Hn being harmonic numbers and ^l>k the polygamma function. One can check 
that hit) = 1 + 0{t^). For the next approximation, i.e., for (^^^"^(a;, t), the nor- 
malization integral is hit) = 1 + Oit*), etc. For ip(^N^oo)ix,t), the normalization 
integral In it) will tend to 1 for all t. 

The normalization integrals versus t are shown in Fig. 2. For approximations 
involving ^oix) and $2(2;), the calculations were done analytically, while the curve 
corresponding to inclusion of ^^™^ix,t) was calculated numerically. As seen from 
Fig. 2, adding more terms brings normalization closer to 1. 



Normalization 
I A 

1 . T-- . -. ; :—. 

0.6 



0.2 

qI ^ J t 

0.5 1 



Fig. 2. Normalization factor calculated for terms including only $0(2;) (dot-dashed line), ^o{^) 
and ^2{x) (dashed line) and ^o{x), ^2{x) and ^:i{x) (solid line). 



In this situation, it makes sense to introduce the "normalized Ansatz" , in which 
fix,t) is approximated by the ratio i'jv(x,t) = Lp(M)ix,t) / I^it), so that the correct 
normalization of the A'^th approximant is guaranteed for all t. In particular, this 
gives 



tTi2 + 2t) l-t^lnxlni 



As seen from this formula (and also from Fig. 3), the initial flat function immediately 
(for whatever small positive t) evolves into a function vanishing at the end points 
with its shape dominated by the (xa;)* factor. 
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Fig. 3. Evolution of the flat DA under the singular part of the evolution kernel: the curves shown 
correspond to t = 0, t = 0.3, t = 0.6, t = 1.0. 



2.3. Adding Non-Singular Part 



When the whole QCD evolution kernel is taken into account, we have 

v{x) = ~ [ [V{y,x) ~V{x,y)]dy ^ 3/2 + xlnx + xlnx (25) 
Jo 

= w™s(a;) - l/2-a;lnx-xln5 . 

Following the same steps as in the previous section, we calculate the expansion 
terms for the initial flat distribution amplitude: 

M^) =1 (26) 

*i(x) -0 (27) 

11 X 

^2{^) —xlnx + x\nx -{ — a; In x -\ — xln x + xx\n Inxlnx 

2 2 X 

+ xLi2 (-|^) +iLi2 (-|) . (28) 
With these terms taken into account we have 

^(2)(x,0=e3*/2 {x'xn' (l + |<i>2(^) 



(29) 



Unfortunately, for this form it is impossible to analytically calculate the normaliza- 
tion integral even for the lowest term. Compared to the Ansatz used for the singular 
part of the evolution kernel, the Ansatz 

ip{x, t) = e*'"(")$(:E, t) = e3*/2(a;^i")*$(a;, t) (30) 
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(a) (b) 
Fig. 4. Expansion components in the full kernel case: ^i{x) (a) and ^2{x) (b)- 



has an extra overall factor [e~^^'^x~^x~^Y . Note that the function x~^x~^ is finite 
at the end points x = 0, 1, where it takes its minimal value for the interval [0, 1] 
(equal to 1), and has a maximum for x = 1/2, where it equals 2. Thus, the factor 
^~i/2^-x^-x gjijiances the xx profile in the middle (by 2/^/e w 1.2 factor) and 
suppresses it at the end points (by l/^e « 0.6). This is a rather mild modification, 
and what is most important, it does not change the ^ x* (or ~ x*) behavior at the 
end points. So, it makes sense to use the expansion 



^(-l)"(a;lnx + SlnS)' 

n=0 



(31) 



in powers of t and combine it with the expansion for <i>(a;,t). This corresponds to 
Ansatz 

ip{x,t) = ixxye''/^ (^^„ix) + t^,{x) + ^^^^2{x) + ..)j , (32) 

whose expansion coefficients $„(a:) can be straightforwardly obtained from $„(a;)'s. 
In particular, <i>i(x) = —{xlnx + xlnx), and $2(2^) ~ ^2{x) + {xlnx + x\nx)^. 

Now, the normalization integral for the lowest terms can be calculated analyti- 
cally: 



t 



t 



m ^ I .(t)(x, t).. . e-/^Llld-| (1 - ^ + ^^i-H. + i^t..) ) ,(33) 

where H„ are harmonic numbers. Fig. 5 shows the normalization versus t. 

Again, we may switch to the normalized Ansatz formed by the ratio 
LP(^[^-^{x,t) / Ij^{t). For a flat initial distribution, this gives 

^V{2 + 2t) l-t{x\nx + x\nx) 

The results are illustrated by Fig. 6. 
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Fig. 5. Same as in Fig. 2 for tlie case of 
the full kernel. 



Fig. 6. Same as in Fig. 3 for the case of 
the full kernel. 



3. Evolution of Anti-Symmetric Flat DA 



3.1. Singular Part 



Evolution equations may be applied also in situations when the distribution ampli- 
tude is antisymmetrie with respect to the change a; — > \ — x. An interesting example 
is the Z^-term function d{x) that appears in generalized parton distributions. Thus, 
let us consider evolution of the DA that initially has the form 



1 < X 1/2 
-1 1/2 < a; < 1 



The v{x) function (111 is the same, since it depends on the kernel only. Thus we 
can use the same Ansatz ([6]) and expansion (|8]). Since ipa{x) is not just a constant, 
the first expansion coefficient $i(a;) is nonzero. Let us start with the singular part 
of the kernel. Then we get 




-l+2x 



< x 1/2 , 
1/2 < x 1 . 



We see that there are logarithmic terms In |1 — 2x| singular for x — 1/2. These terms 
are natural, since each half of the antisymmetric DA on its interval is expected to 
evolve similarly to a flat DA on the < a; < 1 interval. This observation suggests 
the Ansatz 

ip{x,t) = e'^\xxY\l - 2xp*$(a;,t) . (35) 
With this definition of $(a;,t), the ln|l — 2x\ terms are eliminated from $i(a;): 



$i(x) = -21ni 61(0 < X < 1/2) ~ {x ^ x) 



(36) 
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For the expansion component <i>2(a^), we have 



$2(x) =61(0 < a; < l/2)< 



2 In x In a; — In (1 — 2x) — 3 In x(ln(a;a;) — 1) 



21n(l - 2x){\n{4:xx) + \n{xx)) - 8Li2(l - 2x) - 4U2{x) - 4Li2(2a;) 



2Li2 



4Li2 



2x- 1 



4Li2 



1 - 2x 



{x ^ x} . 



(37) 



The graphical resuhs for the expansion components are shown in Fig. 7. 




X) 




(a) (b) 
Fig. 7. Expansion coefficients for the antisymmetric DA; <I>^"'^(x) (a), ^^^^i^) 0^)- 



The evolution of ipsing{x,t) to this accuracy can be obtained from 



ip{x,t) = e^\xxy\l - 2xr ( Mx)+t^iix) + ^*2(a:) 



(38) 



As shown in Fig. 8, the initial step function evolves into a function which is zero 
at the end points and in the middle point. 
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X 



Fig. 8. Evolution of antisymmetric DA under tlie singular part of the kernel for i = 0, f = 0.3, 
t = 0.6, i = 1 



3.2. Adding the N on- Singular Part 

Since the nonsingular part does not add In a; and Inx terms to v{x), we may proceed 
with the same Ansatz 



Lp{x,t) = e3*/2(ra)*|l-2a;|2*$(a;,i), 



(39) 



but the expansion components change (see Fig. 9). 



$i(x) 



1 

0.5 


-0.5 

-1 







0.5 



1 




(a) (b) 
Fig. 9. Expansion components (a) and '^2{x) (b) in the case of the full kernel. 
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The distribution amplitude is now built using 

ip{x, t) = e^/^\xxY\l - 2xp* (^(^o(x) + + ^$2(a;)^ , (40) 

For the first coefficient we have 

$i(a;) 6'(0 < X 1/2) [-2a;ln2 - 21nx - (xlnx + a;lna;)]-{x x] , (41) 
and for the second, 

$2(2;) =61(0 <a;^ l/2)i ^(a:-2) + 21n^ { ^ ^ \ + (3a; - 2) In^ 2 
12 \ 1 — 2x / 

13 \ 

y - 3a; j ln2i + ln(l - 2a:)(4(3 - 2a;) In 2 + 2(5 - 2a;)lnx) 

lnx(l - 2(1 - 3a;) ln2 + (3x - 5) Ina;) + a; ^41na;ln2 + ln 

4xLi2 f j - 2(a; - 4)Li2 f — ^ J - 2(5 + a;)Li2(a;) 

■ 2(5 - 2x)Li2(2a;) - (1 + 2a;)Li2 (-^ 



X 



X \ . X 



+ Li2^^j+4Li2^-— (42) 

As may be seen from Fig. 10, the resulting curves are rather close to those obtained 
when only singular part of the kernel was taken into account. 




X 



Fig. 10. Evolution of the antisymmetric DA in the fuU kernel case. The curves for t = Q, 
t = 0.3, t = 0.6, t = 1 are shown. 
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4. Summary 

We described a new analytical method of performing ERBL evolution. Unlike the 
standard method of expansion in Gegenbauer polynomials, the method works for 
functions that do not vanish at the end points. The method was applied for two 
cases of the initial DA: for a purely flat DA, constant in the whole < x < 1 interval 
and for an antisymmetric DA which is constant in each of its two parts < a; < 1/2 
and 1/2 < a; < 1. In case of a purely flat DA, the leading term gives (a;x)' evolution 
with the change of the evolution parameter t. For the accompanying factor, two 
further terms in the expansion were found. In case of an antisymmetric flat 
DA, there is an extra factor |1 — 2a;p* that takes care of the jump in the middle 
point X — 1/2. The correction terms were also calculated. The results show good 
convergence for t < 1/2. It should be noted that for t > 1/2, the evolved DA is 
rather close to the asymptotic form, and one can use the standard method of the 
Gegenbauer expansion which is well convergent for such functions. 
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